Laplace Transform of Spherical Bessel Functions 

A. Ludu 

Department of Chemistry and Physics, 
Northwestern State University, Natchitoches, LA 71497 

R. F. O'Connell 
Department of Physics and Astronomy, 
Louisiana State University, Baton Rouge, LA 70803-4001 

Abstract 

We provide a simple analytic formula in terms of elementary functions for the 
Laplace transform ji(p) of the spherical Bessel function than that appearing in 
the literature, and we show that any such integral transform is a polynomial of 
order I in the variable p with constant coefficients for the first I — 1 powers, and 
with an inverse tangent function of argument l/p as the coefficient of the power I. 
We apply this formula for the Laplace transform of the memory function related 
to the Langevin equation in a one-dimensional Debye model. 
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I. Introduction. The memory function in the Debye model. 



In recent years, there has been widespread interest in dissipative problems arising 
in a variety of areas in physics, and we refer to [1] for a mini-review. As it turns out, 
solutions of many of these problems are encompassed by a generalization of Langevin's 
equation to encompass quantum, memory and non-Markovian effects, as well as arbi- 
trary temperature and the presence of an external potential V(x). As in [2], we refer 
to this as the generalized quantum Langevin equation (GLE): 

mx + f dt'ii(t-t')x(t') + V'(x) = F(t) + f(t), (1) 

where V'(x) = dV(x)/dx is the negative of the time-independent external force 
and ji(t) is the so-called memory function. F(t) is the random (fluctuation or noise) 
force and f(t) is a c-number external force (due to a gravitational wave, for instance). 
In addition (keeping in mind that measurements of Ax generally involve a variety of 
readout systems involving electrical measurements), it should be strongly emphasized 
that "- the description is more general than the language -" [2] in that x(t) can be a 
generalized displacement operator (so that, for instance, Ax could represent a voltage 
change). Furthermore, fj,(t) and F(t) are given in terms of the parameters of the heat 
bath only. Explicitly 

fi(t) = '^rrijUj cos((Vjt)6(t), (2) 
j 

where 9(t) is the Heaviside step function. Also 

F(t) = £mX<#)> ( 3 ) 

j 

where q h (t) denotes the general solution of the homogeneous equation for the heat- 
bath oscillators (corresponding to no interaction). Thus, we have all the tools necessary 
for the analysis of any heat bath. As emphasized in Refs. [1,3] of primary interest for the 
calculation of observable physical quantities are either the Fourier or Laplace transforms 
of fj,(t). In particular, the blackbody radiation heat bath was investigated in detail in 
Refs. [2,3]. 

Here, we wish to consider the one-dimensional Debye model [4] because it leads to 
a result for the memory function which involves spherical Bessel functions. 

In this model one takes all the heat bath oscillators to have the same mass i.e. 
rrij = m for all j. Thus, taking the maximum allowed heat bath frequency to be w^, we 
see that Eq. (2) reduces to 

fj,(t) = m Uj cosujjt = m I dujD{uj)uj 2 cos ut (4) 
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where [4] 

D{u) = A, (5 ) 

TXV 

is the density-of-states. Also, v is the velocity of sound and L is the length of a 
one-dimensional line such that u — vk where Ak = (2,71 /L) is the interval between the 
allowed values of the wave-vector k. It follows that 

fi(t) = / duuj 2 cos ut = - — LJ 3 L {j (u L t) - 2j 2 {uj L t)} (6) 

7TV JO OTTV 

Thus, in order to calculate jl{p) we will turn to our general result given in the next 
section and related to Laplace transform of spherical Bessel functions. For example, for 
the above case, the following expression for the Laplace is obtained 

fi(p) = —uj L p- p 2 tan' 1 (—) ,. (7) 

TXV \ P J 

from a more general formula. For a detailed discussion of this result we refer to Ref . [5] . 



II. Analytic formula for the Laplace transform. 



In order to calculate eq.(7), as well as the Laplace transform of higher order Bessel 
functions for other applications, we introduce in the following an exact formula, much 
simpler than formulas found in literature, in terms of trigonometric functions and poly- 
nomials. 

The spherical Bessel functions are given by (Ref. 6, p. 965) 

m = (^fji H (t), (8) 

where I is a positive integer and J„ are the Bessel functions of the first kind of real 
argument t. The Laplace transform, defined by 

L\j t (t)] = j t (p) = / me-*dt, (9) 
Jo 

can be calculated by using the relation [Ref. 7, p. 182, II (9)] 

p-" ( p 

(P 2 + 1) 2 
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where P^ are the Legendre functions and we have to fulfill the restrictions: fi + v > — 1 
and p > 0. Also, T is the Gamma function defined [6] as T(x) = e~H x ~ x dt for 

positive values of x. 

In the case of the spherical Bessel functions we have v — I + 1/2 and p = —1/2 
and hence the restriction is I > — 1 which is always fulfilled since I — 0, 1, Hence 

3i(p) = yfr(J + i) Q Vj^(p) 
tt r(< + i) M^fi-pf /i i. 3. i_-p\ 

2T(Z + 3/2) VpJ U+P/ V2' 2 ' 2 ' 2 7' 1 ' 

where p = pj \Jp 2 + 1 and p > so that 1 > p > 0. Here F (a, j3; 7; x) is the Gauss 
hypergeometric series defined by the formula [6] 

F(a, /3; 7; x) = X] 



fc=0 



7 (fc) fc! ' 



where a (fc) = a(a + l)...(a + k - 1); a (0) = I. 

The formula provided by eq.(ll) is somehow difficult to use in the cases I > 2, and 
requires tedious calculations and further manipulations such as integration or other 
transforms. However, in the case / = 0, eq.(ll) reduces to 

Since the Gauss hypergeometric series in eq.(12) can be writen in the form (Ref. 6, p. 
1041 II 13) 

F I ?; p) = S ^l, (13) 

V2' 2' 2' V y/p ' v ; 

we find for the Laplace transform of jo(t) = sint/t the usual expression (Ref. 7, p. 152) 

jo = tan" 1 -. (14) 
p 

Another possibility for expressing the Legendre function for the I = case in a simpler 
way is to use the result ( Ref. 6, p. 1008) 



Pj(p) = hmP;\(p) = cos- 1 ^) j ■ (15) 



Finally, one can use a direct approach, based on the fact that the Fourier representation 
of the spherical Bessel functions are definite integrals over Legendre poynomials 
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The following step is to perform the integration for the Laplace transform, obtaining 
the relation j\{p) = i l+1 Qi(ip), where Qi(p) are the Legendre functions of second type 

However, none of these approaches can provide an analytic expression in terms of 
elementary functions for the / > 2 case. Thus, we are motivated to develop a new 
approach based on recursion relations. By differentiating two times the second order 
spherical Bessel function 

/3 \ sint 2cost 
J2 = U " l ) 1 ^' (16) 

and by using recurrently the formula for the Laplace transform of the derivative of a 
function /(*) (Ref. 7, p. 129) 

f'(p)=pf(p)-f(0), (17) 

where the prime denotes the derivative with respect to p, we can calculate the Laplace 
transform of j^.That is we choose for f(t) = sinat/t which has its Laplace transform 
in tables [6,7]. We differentiate f(t) once and we calculate the Laplace transform of 
f'(t) by using Eq.(17). Then, we differentiate one more time and find again the Laplace 
transform of f"(x) by formula Eq.(17). Since we can express j 2 only in terms of the 
functions f(t) and / (x) we can express its Laplace transform in terms of the Laplace 
transforms of /(£), f"(t). We finally obtain 

f3p 2 l\ ,1 3p 

This procedure gives us a hint for constructing a recursion formula for the Laplace 
transform of the spherical Bessel functions of any order. This will provide a much 
simpler form than that one provided by eq.(ll). This approach needs only the derivative 
recursion formula, eq.(17), and the transform of the zero order function. By starting 
from the recursion formulas for the spherical Bessel functions (Ref. 6, p. 967) 

ji+i(p) = -j[{p) + -3i{p) 
P 

2/ 4- 1 

Ji + i(p) = ^Ji(p)-Ji-i(p), l>h (19) 
P 

and by writing them in a convenient form 

21 + 1 . I 

J'm = -—J, + j^ji-u (20) 

we find that any spherical Bessel function can be expressed as a sum of the derivatives 
of different orders of j 

2k<l i\-1k 

Up) = E (f^M, (2i) 



k=0 
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where are coefficients depending on the labels / and k. In order to calculate these 
coefficients we use again the recursion relation eq.(20) and the formula in eq.(21) and 
we obtain the equation 

2k<l , 2k<l-l ; . -i 2k<l+l 

nil) A (l-2k+i) _ I \- r (l-i) .(i-2fc-i) t + I (j+i) .(i-2fc+i) / 99 n 

2^ u k Jo — 2l + \ ^-^ ~ 21 + 1 ^ ' 

fc=0 fc=0 fc=0 

where the superscript on the Bessel function indicates the order of differentiation. Since 
the Wronskian form of any two derivatives of jo is not identical to zero, all these 
derivatives are functional independent so we can determine the numbers in eq.(22) 
by indentifying the coefficients of the same order of derivative of the Bessel function. 
Consequently we find the binary recursion relation 

= ~ 2 4^CS\ (23) 



Z + l 



Ct 1] = ("ir 1 !^, * = 0, 1, ... and cf = 0. (24) 



which provide all k = coefficients 

[2l_ 

(1 + 1) 



and a ternary recursion relation 

which provides the other coefficients. Also, if we are given Cj 2) (which can be easily 
obtained by giving particular values in eq.(22), that is C[ 2 ^ = 1/2) we can generate all 
the coefficients by using the ternary relation, eq.(25). By using this procedure it is easy 
to identify all coefficients in the formula eq.(21). For example we have 

_ ■(!) ■ _ 3 '(2) , 1 ■ ■ 5 -(3) 3 .(i) 

Ji — Jo > 32 — —Jo ' 2^°' -? 3 ~~ 2^° ' 

k = fj'o 4) + ^jf + \h, etc. (26) 
The next step is to calculate the Laplace transform by using eqs.(17) and (21). Since 

^(0) = ^, J'f +1) (0) = 0, (27) 
and from eqs.(14) and (17), we have 

L \jo n) ] =p"tan -1 - - x!p n_m_1 jo m) (0), (28) 



V 



m=0 
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we finally obtain the desired formula 



Up) =p 1 



1 



2k<l 



tan- - . E c^p 

V k=Q 



(0„-2fc 



2k<l 
fc=0 



(0 



-1^ 



-2fc-2m-l 



(29) 



valid for / > 1, where [ • ] stands for the integer part. The Laplace transform is 
hence given in terms of elementary functions, by a polynomial Q of order I — 1 in 
the variable p, with constant coefficients (all coefficients determined by the recursion 
relations eqs.(24) and (25)) plus another polynomial P of order / multiplied by the 
inverse tangent function of 1/p 



j l = P l (p)t(m- 1 (l/p) + Q l _ 1 (p). 



(30) 



Finally, we compare the above algorithm, based on recurrsion relations Eqs.(24, 25) 
and Eqs. (29, 30), with direct numerical algorithms used in the corresponding built-in 
functions in different symbolic programs, like for instance Mathematica. All functions 
involved above have their equivalent built-in correspondent in Mathematica- 4-0 fused 
on a Power Macintosh 8500 computer), that is Gamma[z], BesselJ[n,z], LegendreP[n,z], 
LegendreQ[n,m,z], etc. While comparing the CPU time elapsed for both procedures we 
find out that our recurrsion approach is faster. That is the relative difference between 
the two CPU intervals of time increases as e 3 ' 21 (numerical with respect to analytic), 
where / is the order of the spherical Bessel function. 



III. Conclusions 



By using a generalization of the quantum Langevin's equation in terms of memory 
function one can obtain exact solutions for dissipative problems arising in many areas in 
physics. In terms of a one-dimensional Debye model the memory function is expressed 
as a combination of spherical Bessel functions. However, of primary interest for the 
claculation of physical quantities are the laplace (or Fourier) transform of the memory 
function. 

Also, the Schrodinger equation for a free particle in polar coordinates leads, for each 
value of the positive integer / of the orbital angular momentum, to a radial equation 
which results in the generic equation for spherical Bessel functions. In electrodynamics, 
too, the spherical Bessel functions are related to solutions of the field equations in the 
stationary or quasi-stationary regime in cylindrical geometry. And moreover, Bessel 
functions are involved in the Helmholtz equation in cylindrical coordinates. In all such 
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applications one needs the integral transforms (especially the Fourier or Laplace trans- 
forms) of such solutions. Constructing an exact analytic formula for such applications 
is the object of our studyy. Consequently, the present letter is not addressed to experts 
in special functions, but rather to physicists who just need to apply such formulas in 
their research, in the simplest possible form. 

We have obtained a simpler exact formula for the Laplace transform of the Spherical 
Bessel functions of any order, in terms of polynomial and trigonometric functions, that 
is in terms of elementary functions. Among other representations of such functions, this 
expression is much simpler, hence more useful for potential applications. The speed of 
calculation is one of foremost importance for applications similar with the one presented 
above, especially in more than one dimension. We compared this formula with some 
numerical evaluation in terms of precision and suitability for faster implementation. 
Our comparison shows that the introduced formula is easy to be used and it introduces 
less computational complexity than conventional numerical techniques. 
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